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ABSTRACT 


Analytical  studies  of  the  hole-pressure  error  for  non-Newtonian  creeping  flows  over  a  trans¬ 
verse  slot  are  pursued  with  particular  interest  in  the  theory  of  Higashitani,  Pritchard,  Baird  and 
Lodge  (HPBL).  To  correctly  apply  the  HPBL  theory  a  modified  hole-pressure  relation  (MHPR) 
is  employed.  Some  important  mathematical  properties  of  the  MHPR  are  presented.  By  studying 
the  MHPR  in  streamline  coordinate  formulation  we  find  a  fortuitous  error  cancellation  phe- 
nomenon  in  the  derivation  of  HPBL  formula.  For  second-order  fluid  and  Tanner’s  viscometric 
moder  (under  certain  assumptions)  the  error  cancellation  is  proved  to  be  exact.  It  is  this  er¬ 
ror  cancellation  phenomenon  that  provides  a  complete  theoretical  explanation  for  the  paradox 
between  an  apparently  flawed  derivation  and  the  fortunate  success  of  the  HPBL  prediction. 
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1.  Introduction 

In  this  paper  we  consider  shear  flow  of  some  liquids  between  parallel  walls  which 
axe  at  a  separation  of  h.  One  wall  contains  a  hole  whose  cross  section  is  a  transverse 
slot  (of  narrow  dimension  w  parallel  to  the  main  flow  direction).  There  is  no  net  flow 
through  the  hole.  A  transducer  at  point  a  on  the  bottom  of  the  hole  records  a  normal 
thrust  Pa  =  — T22U,  exerted  by  the  fluid.  Another  transducer  mounted  flush  with  the  wall 
face  opposite  the  hole  records  a  normal  thrust  Pj,  =  —  T22I&  at  a  location  b  on  the  hole 
centerline  (£  shown  schematically  in  Figure  1,  is  a  well-established  experimental  procedure 
for  determining  the  fluid  thrust  component  normal  to  a  solid  wall  in  the  undisturbed  shear 
flow.  For  slow  flows  of  Newtonian  fluids,  it  is  well  known  that  the  two  readings,  Pa  and 
P&,  sire  the  same,  i.e.  no  error  is  introduced  by  this  procedure.  But  for  non-Newtonian 
fluids,  by  contrast,  the  difference  between  Pa  and  Pi  can  be  very  substantial  as  shown  both 
exp erimet ally  and  numerically  (see,  among  others,  references  [1]  to  [12]).  This  difference 
is  called  hole-pressure  error  in  the  literature  and  usually  denoted  as  P//,  namely 

PH  =  Pi  -  Pa  =  -Of22!i  -  f22|a)  .  (1.1) 

The  physical  explanation  of  the  existence  of  hole-pressure  error  Ph  can  be  found  in  refer¬ 
ences  [13,14].  Roughly  speaking,  the  stretching  of  the  N\-  spring,  as  shown  in  Figure  1, 
indicates  how  an  extra  tension  T  along  a  streamline  lifts  the  fluid  out  of  the  hole  resulting 
in  a  low  reading  of  Pi-  An  analytical  calculation  of  the  hole- pressure  error,  Ph,  was  first 
made  by  Tanner  and  Pipkin  [2]  for  the  particular  case  of  slow,  two-dimensional  flows  of 
second-order  fluids  past  a  transverse  slot.  They  showed  that  Ph  is  equal  in  magnitude  to 
a  quarter  of  the  first  normal  stress  difference,  Ni,  of  the  fluid. 

The  hole-pressure  problem  has  been  a  research  topic  of  considerable  interest  in  recent 
years  because  of  the  possibility  of  obtaining  elasticity  measurements  for  viscoelastic  liquids. 
Of  particular  interest  has  been  the  theory  of  Higashitani,  Pritchard,  Baird  and  Lodge 
(HPBL)  [15,5,6]  which  leads  to  expressions  relating  Ph  to  the  shear  stress  Ti2,  the  wall 
shear  stress  aw  and  the  first  normal  stress  difference  N j.  Based  on  the  phenomenon  of 
hole-pressure  error,  Lodge  and  de  Vargas  [5]  constructed  a  slit  die  rheometer  using  the 
hole-pressure  difference  P#  to  measure  the  first  normal  stress  difference  N 1  for  molten 
polymers.  This  application  has  made  it  more  important  to  understand  the  theory  behind 
the  hole-pressure  phenomenon,  and  the  hole-pressure  problem  itself  has  been  proved  to  be 
a  highly  controversial  issue. 

The  big  conflict,  that  has  been  puzzling  people  in  this  field  for  many  years,  is  that 
although  very  good  agreement  between  the  theory  of  HPBL  and  the  published  experimental 
&  numerical  results  has  been  obtained  for  several  polymer  solutions,  polymer  melts  and 
some  non-Newtonian  constitutive  models,  it  has  been  known  that  at  least  two  of  the 
key  assumptions  in  the  theory  are  violated  [16].  Some  attempts  have  been  made  in  the 
literature  trying  to  correct  the  flaws  in  HPBL  theory  and  to  explain  the  paradox,  but  not 
much  progress  has  been  achieved  so  far. 
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Figure  L  Schematic  illustration  of  hole-pressure  prob¬ 
lem  and  definition  of  streamline  coordinate  system  for  pla¬ 
nar  flows  past  a  transverse  slot. 


In  order  to  correctly  apply  the  HPBL  theory,  a  modified  hole-pressure  relation 
(MHPR)  in  path  integral  form  is  used  in  our  work.  This  MHPR  contains  an  extra  term,  i.e. 
the  contribution  of  the  stress  gradient  This  term  is  non-zero  in  general,  according 

to  our  numerical  observation  on  the  excess  pressure  rise  phenomenon  [17],  however  it 
was  assumed  to  be  zero  and  was  neglected  in  the  original  derivation  of  HPBL  theory. 
To  investigate  the  two  key  flaws  of  the  HPBL  theory,  we  found  that  it  is  sufficient  for 
the  present  to  consider  the  second-order  fluid,  since  both  flaws  exist  there.  By  studying 
the  MHPR  based  on  the  streamline  coordinate  formulation,  we  found  a  very  fortuitous 
and  interesting  phenomenon  behind  the  success  of  HPBL  prediction  :  the  exact  error 
cancellation  for  the  second-order  fluids ,  i.e.  the  error  caused  by  one  key  flaw  can  be 
exactly  canceled  out  by  the  error  introduced  by  another  key  flaw.  By  using  the  streamline 
coordinate  formulation,  we  are  also  able  to  prove  analytically  the  equivalence  of  the  MHPR 
with  the  re-interpreted  HPBL  prediction.  That  is,  although  the  two  relations  have  different 
forms  and  contain  different  terms,  the  MHPR  and  the  re-interpreted  HPBL  theory  predict 
the  same  exact  value  of  N®w/ 4  for  the  second-order  fluid.  Furthermore,  the  conclusion 
of  error-cancellation  has  also  been  generalized  to  include  Tanner’s  “viscometric  model.” 
These  theoretical  results  enable  us  to  eventually  solve  the  puzzle  about  the  success  of 
HPBL  prediction  and  to  explain  the  controversy  in  the  derivation  of  the  HPBL  formula. 

The  outline  of  this  paper  is  as  follows.  In  section  2  we  review  the  HPBL  theory, 
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some  recent  work  in  the  literature  and  the  modified  hole-pressure  relation;  in  section  3  we 
develop  the  streamline  coordinate  formulation  of  the  MHPR  for  a  second-order  fluid;  in 
section  4  we  prove  the  exact  error  cancellation  in  MHPR  for  second-order  fluid  and  for 
Tanner’s  viscometric  model;  finally,  in  section  5  we  discuss  our  conclusions. 

2.  The  HPBL  Theory  and  Modified  Hole-Pressure  Relation 
2.1  HPBL  Theory  of  Creeping  Flow 

By  considering  the  equation  of  motion  in  streamline  coordinate  form  and  making 
some  assumptions  about  the  flow  field,  Higashitani  and  Pritchard  [15]  deduced  a  proposed 
relation  between  hole-pressure  difference,  P #,  and  the  first  normal-stress  difference,  N\. 
The  final  relation  is  attributable  to  Higashitani,  Pritchard,  Baird,  and  Lodge  [5,6],  and  it 
has  two  forms  for  creeping  flow: 

academic  form 


practical  form 


Ph  = 


? 


Nlw  =  2  PH 


din  Ph 
diner  w 


or  Niw  =  2  at 


dPH 

da 


(2.1) 


(2.2) 


where  aw  and  N\w  are  the  values  of  T\2  and  N\  at  the  channel  walls  as  shown  in  Figure 
1.  Ph  is  the  hole-pressure  difference  defined  by  equation  (1.1).  The  practical  form  (2.2) 
is  obtained  by  differentiating  (2.1)  with  respect  to  aw.  It  is  important  to  note  that  in  the 
original  HPBL  formula  N\  &  r  correspond  to  the  steady  unidirectional  shear  flows  only, 
and  N\  =  v\i2  &  r  =  777  are  the  strict  viscometric  functions.  However,  in  this  paper 
we  shall  define  N\  as  the  generalized  viscometric  part  of  Tn  —  T22  for  more  general  flows 
which  may  not  necessarily  be  a  unidirectional  shear  (viscometric)  flow. 

For  convenience  of  further  discussion  and  reference,  the  assumptions  made  either 
explicitly  or  implicitly  by  Higashitani,  Pritchard,  Baird  and  Lodge  in  their  derivation  of 
equations  (2.1)  &  (2.2)  [15,5,6]  are  summarized  as  following: 

(Al)  Incompressible,  isothermal,  creeping  Couette  flow  (plate- driven)  is  primarily  con¬ 
cerned.  (Re  is  very  small) 

(A2)  The  flow  patterns  axe  symmetric  about  the  hole  centerline  (£.  Particularly,  the 
streamlines  are  symmetric  about  (£,  which  gives  the  symmetry  conditions 

p-=0.  (.  =  1,2)  V(„,92)6(E  (2.3) 

Oqi 

(A3)  An  orthogonal  curvilinear  coordinate  system  based  on  the  streamlines  of  the  flow 
field  can  be  established  in  the  local  area  around  the  hole  centerline. 
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(A4)  The  stresses  axe  symmetrically  distributed  about  the  hole  centerline  (£,  consequently 


dT\ 2 

dqx 


=  0 


V  ( 9i>?2 )  £  € 


(2.4) 


?h=0  V(„,g2)€®  (2.5) 

Oqi 

( A5)  The  motion  at  the  hole  centerline  is  that  of  a  shear  (viscometric)  flow.  Therefore,  the 
quantities  N\  and  Ti2  in  equation  (2.1)  are  the  same  as  the  viscometric  functions. 
(A6)  The  path  integral 


exists  and  the  change  of  variable  from  streamline  coordinate  q2  to  shear  stress  a 
should  be  valid. 

It  is  known  from  [17]  that  the  streamwise  gradients,  and  are  non-zero  on 

<£  even  for  the  Stokes  solution  of  Couette  base  flow  over  a  slot.  Flow  visualization  [6] 
and  analytical  studies  [17]  have  shown  that  Tn  —  T22  is  not  viscometric  on  <£  due  to  the 
existence  of  slot.  These  two  important  facts  suggest  that  the  derivation  of  HPBL  prediction 
(2.1),  especially  when  it  is  desired  to  apply  HPBL  theory  to  Poiseuille  base  flow,  is  flawed 
mainly  in  two  ways: 

(i)  The  assumption  of  Tn,i  =  0  on  <£  is  incorrect,  i.e.  the  streamwise  stress  gradient 
is  in  general  non-zero.  (We  will  see  later  on  in  section  4  that  the  term  Tii,i  plays 
an  important  role  in  the  error  cancellation.) 

(ii)  The  assumption  (A5)  is  not  true,  i.e.  T\\  —  T22  is  in  general  not  viscometric.  There¬ 
fore,  the  change  of  variable  required  in  getting  the  rheological  space  integral  (2.1) 
from  the  path  integral  (2.6)  is  invalid.  The  change  of  variable  requires  that  Tn  —  T2 2 
can  be  written  as  a  unique  function  of  r  =  Ti2(<?2).  This  is  not  possible  in  general 
as  will  be  seen  in  equation  (3.17).  It  is  also  easy  to  see  that  the  correspondence 
between  the  values  of  Tn  —  T22  and  the  values  of  T12  is  not  one  to  one. 

These  two  flaws  have  been  pointed  out  by  Malkus  &  Yao  in  reference  [16].  From  assump¬ 
tion  ( Al)  we  know  that  only  creeping  Couette  flow  is  primarily  concerned  when  Higashitani 
and  Pritchard  derived  the  relation  (2.1)  in  [15].  Since  that  time,  the  possibility  of  applying 
relation  (2.1)  &  (2.2)  to  Poiseuille  flow,  and  to  flows  with  inertial  effects,  etc.  has  been 
explored  both  experimentally  and  numerically  in  references  [5-12,18].  There  is  now  con¬ 
siderable  experimental  and  computational  evidence  to  suggest  that  (2.1)  is  a  reasonably 
good  approximation  to  the  observations. 

Srinivasan,  in  a  recent  paper  [19],  attempted  to  extend  the  HPBL  theory  by  including 
the  inertial  effects  and  the  effect  of  asymmetry  in  flow  field,  and  the  integral  was  theoreti¬ 
cally  evaluated  along  a  generalized  path  of  zero  slope  of  the  streamlines  (this  path  may  not 
coincide  with  the  hole  centerline  in  the  most  general  situation).  Unfortunately,  although 
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Srinivasan  recognized  that  the  assumption  of  (2.5)  would  be  a  significant  source  of  error, 
he  generally  accepted  the  assumption  -g ^  =  0  in  his  formulation. 

2.2  Tanner’s  Recent  Work  -  an  Alternative  Approach 

To  avoid  the  paradoxes  and  problems  with  the  original  HPBL  analysis,  Tanner  used 
an  alternative  approach  in  reference  [20]  to  formulate  the  hole-pressure  difference.  In  his 
work,  Tanner  made  the  following  fundamental  assumptions: 

(Tl)  Assume  the  fluid  is  desrcibed  by  the  “ viscometric ”  constitutive  equation  [21] 

T  =  -pi  +  r)A  +  (vi  +  V2)A.2  -  (2.7) 

where  r/( 7),  17(7),  ^2(7)  are  functions  of  7. 

(T2)  Assume  that  the  normal  stress  coefficient  17  is  proportional  to  the  viscosity  *7(7), 
viz. 

17  =  arj  (2.8) 

where  a  is  a  constant. 

(T3)  The  following  identity 

V  •  h(B  -  A1)]  +  V  ■  =  V  ^  (fij  (2.9) 

is  used  with  the  postulation  that  Dui/Dt  is  small  and  negligible. 

With  the  above  three  assumptions,  a  closed-form  solution  of  the  modified  pressure  can  be 
obtained,  viz. 

P  =  P°  ~  +  \  2/17^7^  +  z'272  (2.10) 

where  p°  is  the  pressure  of  the  generalized  Newtonian  flow,  i.e.  the  flow  that  has  the  same 
symmetry  as  Stokes  flow  but  its  viscosity,  77(7)  can  be  more  general.  By  calculating  the 
hole-pressure  difference  directly  via.  (2.10),  Tanner  found  the  formula 


(2.11) 


When  the  flow  is  in  the  second-order  regime,  equation  (2.11)  gives 

[2] 

Pti  -  N\w/4  ■ 


the  well  known  result 
(2.12) 


Based  on  this  result,  Tanner  concluded  that  the  HPBL  equation  is  applicable  to  the  hole- 
pressure  problem,  even  though  when  flow  is  not  strictly  a  viscometric  flow  assumed  in 
(Tl).  Furthermore,  Tanner  predicted  that  the  Maxwell  and  Oldroyd-B  models,  which  have 


6 


constant  viscosity  and  normal  stress  coefficients  in  truly  viscometric  flows,  are  expected  to 
show  the  second-order  fluid  result  (2.12). 

Tanner’s  work  is  important  in  verifying  and  supporting  the  conclusion  about  the 
validity  of  the  HPBL  theory  from  an  alternative  approach.  But  his  work  did  not  directly 
resolve  the  paradoxes  and  the  problems  in  the  derivation  of  HPBL  theory. 


2.3  Modified  Hole-Pressure  Relation  (MHPR) 

If  one  goes  through  the  details  in  deriving  equation  (2.1)  by  following  the  way  described 
in  [15],  one  would  find  that  the  equation  in  streamline  coordinates  which  is  integrated  along 
the  hole  centerline  to  give  the  thrust  difference,  Ph,  is  not 

df12  2Ti2  df22 
dq2  fn  -  T22  dq2 


=  0 


(2.13) 


but 


1  dTu  1  dT\2  2  Tj2  <97*22 
hi  dqi  h2  dq2  h2fn-  t22  dq2 


(2.14) 


where  T)j  are  the  physical  components  of  stress  in  streamline  coordinates  and  ht  are  the 
metric  coefficients.  The  first  term  in  (2.14)  is  the  contribution  of  the  streamwise  stress 
gradient.  It  is  assumed  to  be  zero  in  the  HPBL  theory  by  assumption  (A4).  If  we  leave 
the  stress  gradient  term  in  (2.14)  and  continue  the  derivation,  we  will  get  a  modified 
hole-pressure  relation  [16] 


Ph 


fq *  /  1  dTn  1  dfl2\ 
qa  y  hi  dqi  h2  dq2  J 


Tji  ~  T22 

2T12 


h2  dq2  , 


(2.15) 


which  looks  quit  different  from  the  HPBL  predition  (2.1).  (2.15)  is  in  path  integral  form 
and  contains  an  extra  term,  while  (2.1)  is  an  integral  in  the  rheological  stress  space.  Since  it 
has  been  known  that  the  HPBL  prediction  (2.1)  is  quite  successful  for  some  non-Newtonian 
fluids,  and  it  also  seems  to  work  even  under  less  restrictive  assumptions  [5,6,17,18,24,25], 
the  questions  which  arise  are:  What  is  the  relation  between  MHPR  and  HPBL  predic¬ 
tion?  How  does  the  MHPR  (2.15)  work?  Does  (2.15)  predict  the  same  result  as  the 
HPBL  theory?  We  shall  try  to  answer  these  questions  in  the  following  sections. 


2.4  Some  Mathematical  Properties  of  MHPR 

We  have  skipped  over  some  important  details  in  getting  the  path  integral  (2.15)  from 
(2.14).  We  need  to  look  into  the  derivation  more  deeply  and  try  to  define  the  MHPR 
rigorously. 

After  considering  the  symmetric  conditions  (2.3)  &  (2.4),  the  equations  of  motion  axe 
simplified  as  follows 

1  dTn  1  df\2  2T’12  dhi  _  ,9  .gx 

h\  dq\  h2  dq2  hih2  dq2 
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Figure  2.  An  example  of  the  possible  change  of  shear  stress  T12  on 
the  hole  centerline  (£  for  a  Poiseuille  base  flow  of  second-order  fluid  over 
a  slot  (De  «  0.2).  Note  q2  =  0  is  the  center  of  channel. 


1  dh\  1  dT22 

h\  dq2  Ti\  —  T22  ^?2 


(2.17) 


To  obtain  (2.15),  we  eliminated  dhi/dq2  from  (2.16)  &  (2.17),  multiplied  both  sides  of 
(2.14)  by  h2(Tn  —  T22)/2Ti2,  and  then  integrated  it  on  <£.  By  doing  this,  we  have  intro¬ 
duced  singularities  to  the  path  integral  of  MHPR,  because  T\2  can  be  zero  at  some  places 
on  <£.  Figure  2  shows  an  example  of  the  possible  change  of  Ti2  on  <£  for  a  Poiseuille  flow 
of  second-order  fluid.  We  can  classify  the  places  where  Ti2  =  0  into  two  different  cases:  a 
closed  interval  such  as  the  part  [q%,  q%  ]  shown  in  Figure  2;  or  some  isolated  points,  like 
the  points  q2' ,  ( i  =  1,2)  in  Figure  2.  These  two  cases  need  to  be  considered  separately 
and  the  results  are  given  in  the  following  two  propositions. 

Proposition  2.1  For  incompressible,  steady,  creeping,  2-D  flows  over  a  slot  (Re  =  0), 
if  the  shear  stress  Ti2  =  0  on  a  closed  interval  [q2,q2  ]  of  £>  then  the  hole-pressure  difference 
between  q%  and  q%  is  zero,  i.e. 


/  1  dtu 

\hi  dqi 


1  dfi2\ 
h2  dq2  J 


Tn  —  f22 

2T12 


h2  dq2  =  0  . 


(2.18) 


Proof :  Since  Ti2  =  0  on  [q2,q2],  we  have  =  0  on  [q$,  q$]  according 
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to  (2.16).  This  immediately  leads  to  the  result  in  equation  (2.18).  0 

Proposition  2.2  Let  q%  be  an  isolated  point  on  <£  where  T12  =  0.  Assume  fn  and  Ti2 
are  C2  at  point  q Then  the  integrand  of  (2.15)  has  only  a  removable  singularity  at  and  its 
limit  at  q\  exists. 

Proof :  T12  =  0  at  q2  implies  that 


1  dfn  1  dfn 

hi  dqi  /i2  dq2 


at  q2 


from  equation  (2.16). 


Do  the  Taylor  series  expansion  in  the  vicinity  of  point  q% ,  we  obtain 


1  dfn  1  dfi2 

h  i  dq\  h2  dq2 


_d_ 

dq2 


(  1  dfn 

\hi  dqi 


1  dT\2 

h2  dq2 


Aq2  +  0(|A92|2) 

9j 


fi2  = 


dfn 

dq2 


|g.  A q2  +  0(|Ag2|2) 


where  A  g2  =  <?2  ~  • 

gives 


lim 

?2— 


Substitution  of  the  above  two  equations  into  the  integrand  of  (2.15) 


1  dfn  1  df n 


(2.19) 


(Tii  —  T22),j  <  M 


where  constant  M  is  a  finite  real  number.  0 


Remark-.  Since  the  singularity  of  the  integrand  of  (2.15)  is  removable  at  the  isolated 
pole  <72,  the  integrand  will  become  a  continuous  function  at  q £  if  we  define  its  value  to  be 
its  limit  (2.19).  In  this  way  there  will  be  no  real  singularity  involved  in  the  MHPR  and 
the  path  integral  (2.15)  can  be  well  defined.  Based  on  Propositions  2.1  and  2.2,  one  can 
directly  prove  the  following  proposition  regarding  the  existence  of  MHPR  path  integral. 

Proposition  2.3  Assume  Tn  and  T\2  be  C2  continuous.  Assume  fi2  =  0  on  [$2,g|  ] 
and  at  some  isolated  points  q (i  =  1,2,.. .).  If  we  interprete  the  path  integral  of  MHPR  (2.15) 
as 


Ph 


(  1  dfn 
\hi  dqi 


1  ^T*i2 

dq2 


Tn  —  T2  2 
2fi2 


/i2  dq2 


(2.15a) 


where  the  values  of  integrand  at  the  removable  poles  q% ,  qf,  (i  =  1, 2, . . .),  should  be  replaced 
by  its  limits  (2.19),  then  the  path  integral  (2.15)  exists  and  is  finite. 

Remark:  In  practice,  it  may  not  be  necessary  to  assume  as  much  smoothness  as  we 
have  to  obtain  the  same  results.  For  example,  the  C2  continuous  condition  for  fn  &  fn 
assumed  here  may  be  replaced  by  some  weaker  assumption  if  the  distribution  theory  is 
used. 
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Finally,  as  pointed  out  by  Higashitani  &  Pritchard  in  [15],  the  path  integral  in  (2.15) 
may  pass  through  some  points,  such  as  the  center  of  a  vortex,  at  which  the  streamline 
coordinate  system  is  not  well  defined.  The  path  integral  should  then  be  carried  out  in  a 
piecewise  manner. 

3.  Formulation  of  the  MHPR  for  a  Second-Order  Fluid 

In  this  section  we  shall  confine  our  attention  to  the  slow  flows  of  a  second-order  fluid 
and  try  to  formulate  analytically  the  MHPR  in  terms  of  the  Stokes  solution.  First  it  is 
worth  noting  that  the  MHPR  path  integral  in  equation  (2.15)  is  based  on  the  streamline 
coordinate  system;  the  stresses  involved  are  the  physical  stress  components  in  streamline 
coordinates  and  d/dq\  is  the  streamwise  derivative.  An  important  fact  we  will  use  in  our 
development  is  that  for  the  assumed  symmetric  velocity  field  the  streamline  coordinate 
systems  on  <£  are  orthogonal  curvilinear  systems  (by  assumptions  (A2),  (A3)  in  section  2 
and  by  choosing  the  <?2  axis  properly),  and  these  systems  coincide  with  the  local  rectangular 
Cartesian  coordinate  systems.  As  defined  in  reference  [22],  the  physical  components  of  a 
vector  or  second-order  tensor  at  a  point  P  relative  to  a  system  of  orthogonal  curvilinear 
coordinates  are  simply  the  Cartesian  components  in  a  local  set  of  Cartesian  axes  tangent  to 
the  coordinate  curves  through  P.  Any  tensor  equation  involving  only  tensor  values  at  one 
point  and  no  derivatives  with  respect  to  the  coordinates  has  the  same  form  in  terms  of  the 
physical  components  in  orthogonal  curvilinear  coordinates  as  it  has  in  terms  of  rectangular 
Cartesians.  However  the  constitutive  equation  we  are  going  to  consider  here  does  involve 
derivatives  with  respect  to  the  streamline  coordinates.  Accordingly,  we  need  to  work  on 
differentiations  with  respect  to  orthogonal  curvilinear  coordinates. 

3.1  Formulation  of  Stress 

By  the  Tanner-Giesekus  theorem,  the  Stokes  solution  for  the  velocity,  u,  with  a  mod¬ 
ified  pressure,  p,  satisfies  the  second-order  fluid  equations  [14,15,23].  The  resulting  stress 
field  is 


T  =  —  [  P°  ~  +  u2p%)}  1  +  rjA 

—  *?r(B  —  A2  —  I72 1)  +  r/Tm{X2  —  72I)  (3.1) 

where  A  and  B  are  the  first  two  Rivlin-Ericksen  tensors;  77,  T  and  T*  are  the  constants 
of  the  second-order  fluid  expansion,  u  and  p°  are  the  velocity  and  pressure  fields  of  Stokes 
flow,  and  j2  =  trA2/2. 

Define  the  streamline  coordinate  system  as  shown  in  Figure  1.  Particularly,  we  can 
choose 

qi  =  qi(x\,x2),  q2  =  x  2  (3.2) 

^1=^1(91,92),  ^2  =  1-  (3.3) 

From  the  symmetry  assumption  of  the  flow  field,  it  is  not  difficult  to  prove  the  following 
symmetry  conditions: 
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(a)  The  (jq-component  of  velocity,  tii(?i,52)>  Is  an  even  function  of  qi,  consequently 

1  du\ 


,  ~  =0  on  <£ 

hi  dq i 

(b)  There  is  no  flow  across  the  streamline,  i.e. 

U2  =  0  V  (qi,q2)  €  D, 

(c)  p°(qi,q2)  is  an  odd  function  of  qi,  viz. 

p°  =  0  on  (£ 

(d)  2i2(<7i,<72)  is  an  even  function  of  q\ 

1  <9Ti2 


hi  dqi 

(e)  The  derivatives  of  scale  factors  are 


0  on  <£ 


(3.4) 


(3.5) 


(3.6) 


(3.7) 


|^=0,  =  unknown  #  0 

oq i  oq2 

^1=0  ^  =  0 
dq,  '  a?2 


(3.8a) 


(3.86) 


By  applying  equations  (3.4),  (3.5)  and  (3.8),  we  have  the  strain-rate  tensor  on  <£ 


A  =  2D  = 


0  7 

7  0 


(3.9) 


where  7  is  the  shear  strain  rate  and  defined  by 

1  du\  u  1  dhi 

h2  dq2  h\h,2  dq% 

Hence  we  have 


A2  =  72I  . 

The  material  derivative  of  A  is  reduced  to 

1  d 2  Ui  1  1 _ dh  t 


DA 

Dt 


=  2ui 


-  -  J  |  1  on\  • 

TJ  dqf  t  hxk7  dq2^ 

0 


_  (  -I _ 1 -  dAl. 


(3.10) 

(3.11) 

(3.12) 


In  getting  (3.12),  the  following  relation  is  used 

u  1  d2h2  _  1  52ui 

hjh2  d<?2  h2  dqj 


on  <£ 


(3.13) 
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which  can  be  verified  by  differentiating  the  continuity  equation.  The  convective  term  is 


AVu  = 


0  7 

7  0 

- 

«i 

a 

1 At . 

hi  hi 

dqi 

0 

Ul  dh 


h\  h.2 

y  0 


1  8tti  n 

h  i  dq? 

0 


0 


hi  J 


(3.14) 


Using  (3.12)-(3.15)  we  obtain  the  second-order  Rivlin-Ericksen  tensor  on  (£ 

0 


B  =2 


u\  d  ui 
0 


h\  dq\  +1 


(3.15) 


Inserting  (3.9)-(3.11),  (3.15)  into  (3.1)  yields  the  physical  components  of  Cauchy  stress 
tensor  on  (£ 


T=-^°-TrM]1+ri 


0  7 

7  0 


ui  52«i 

0 

“1  1  -v2 

2  qT 

h\  dql 

0 

+  bT 

72  0 

0  72 

h\  dq*  ^  J 

(3.16) 


From  equation  (3.16)  we  can  obtain  some  important  results: 

I.  On  (£,  the  physical  stress  components 

rp  _  f  1  dui  Ui  dhi\  _  . 

12  ^  \  h2  dq2  h\ h.2  dq2  )  ^ 

(  9  Ui  d2u\\  Ui  32Ui 

Ti,-T22=2,T(72-2ri^i)=lV1-4,T-1  1 


H.  The  hole-pressure  error 


'h2  dq\  J  1  h\  dq\  ) 


(3,17) 


Ph  =  -f22|,.  =  \_r,T(  7°„)2  =  N°J4 


(3,18) 


where  7°  and  N®w  are  the  values  of  disturbed  7  and  N 1  at  wall  &;  (£  opposte  to  the  slot. 
It  is  necessary  to  make  the  following  comments  before  we  proceed  further. 

(1)  The  numerical  results  show  that  the  streamwise  derivative,  d2u\/dq2,  is  non- zero 
on  (£.  Therefore,  the  last  equation  in  (3.17)  indicates  that  the  flow  is  not  strictly 
viscometric,  owing  to  the  disturbance  induced  by  the  slot.  Consequently  the  as¬ 
sumption  (A5)  in  the  original  HPBL  theory  is  inappropriate. 

(2)  For  second-order  fluid  under  the  symmetry  conditions,  the  shear  stress  on  (£  has 
the  same  fo:m  as  that  in  viscometric  flow,  the  flow  is  locally  shear  flow  and  7 
is  the  usual  shear  strain  rate  on  <£.  However  the  shear  flow  on  <£  may  not  be 
unidirectional. 
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(3)  The  notation  N\  in  the  last  equation  of  (3.17)  has  different  meaning  from  convention 
and  here  N\  ^  Tn  —  T22  unless  the  flow  is  strictly  viscometric.  In  this  paper  we 
define  the  first  part  in  the  last  equation  of  (3.17) 


JVi  =  ^i72 


(3.19) 


as  the  generalized  viscometric  part  of  Tn  -  T2 2-  We  call  it  the  uviscometric 
part ”  because  it  has  the  same  form  as  the  conventional  viscometric  function  in 
the  unidirectional  shear  flow.  Note  that  even  when  flow  is  symmetric,  Tn  —  T22 
still  may  not  be  the  same  as  N 1  in  the  viscometric  base  flow  because  in  general 
d2u\/dq\  7^  0  on  <£ . 

(4)  Result  I  applies  when  the  depth  of  pressure-hole,  d ,  is  deep  enough.  This  will  be 
the  case  for  d  =  00  and  approximately  so  for  d  w. 

(5)  It  is  known  that  in  general  N°w  N\w  for  a  given  flow  field.  Strictly  speaking, 
for  second-order  fluid  the  correct  hole-pressure  error  should  be  IV °w/ 4  as  given  by 
(3.18)  and  the  HP  result  of  (2.12)  is  only  an  approximation.  However  the  results 
(2.12)  and  (3.18)  give  only  a  slight  difference  in  the  predicted  values  of  Ph  for  a 
given  flow  field  (see  item  (8),  below). 

(6)  In  order  to  make  the  HPBL  prediction  be  consistent  with  the  directly-calculated 
values  of  Ph,  we  define  the  re-interpreted  HPBL  formula  as 


(2.1') 


where  cr^,  is  the  disturbed  wall  shear  stress  at  point  q\  shown  in  Figure  1.  In  (2.1') 
the  original  HPBL  formula  (2.1)  is  re-interpreted  in  two  places:  first  aw  is  changed 
to  cr^;  second  N\  represents  the  generalized  viscometric  part  of  Tu  —  T2 2-  Note 
that  the  original  HPBL  prediction  (2.1)  is  actually  an  approximation  for  the  re¬ 
interpreted  formula  (2.1').  (2.1')  is  precise  in  theory  and  (2.1)  is  convenient  for 
application. 

(7)  The  slight  difference  between  (2.1')  and  (2.1)  does  not  affect  the  application  of 
HPBL  theory.  (2.1)  and  (2.1')  lead  to  the  same  differential-form  N\-Ph  relation. 
As  a  matter  of  fact,  differentiating  bothsides  of  (2.1')  with  respect  to  cr“  yields 


N 


0 

lu; 


=  2(7 


0 

W 


dPH(<r°J 

daaw 


(2.2') 


here  we  use  Ph(<*w)  tj°  emphasize  that  Ph  is  a  function  of  u®  .  It  is  easy  to  see,  by 
comparing  (2.2')  with  (2.2),  that  once  we  obtain  the  differential  relation  between  N 1 
&  Ph,  either  in  the  form  of  (2.2)  or  (2.2'),  and  apply  the  relation  in  the  rheological 
space,  the  values  of  aw  &  <7^  in  (2.2)  &  (2.2')  become  unimportant,  they  only 
correspond  to  two  neighboring  points  on  the  Ti2-axis  in  the  rheological  stress  space 
and  hence  they  can  be  simply  replaced  by  a  dummy  variable,  say,  <7  or  r.  Therefore 
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(2.2)  and  (2.2')  actually  represent  the  same  Ni-Ph  relation  in  the  rheological  stress 
space  and  we  can  simply  rewrite  them  as 

•  (2.2") 

Obviously,  (2.2")  is  the  more  general  differential-form  expression  for  the  Ni-Pfj 
relation  in  the  rheological  stress  space,  and  (2.2)  &  (2.2')  are  two  special  cases  for 
(2.2")  when  a  takes  the  values  of  crw  ic  ,  respectively. 

(8)  (2.2")  is  only  useful  for  theoretical  analysis.  In  practice  the  values  of  (or  7») 
axe  not  available  experimentally  and  some  approximation  has  to  be  made.  Usually 
the  following  combination  of  (2.2)  and  (2.2')  is  employed,  i.e. 


Nlw 


—  2aw 


dPn(al) 


(2.2'") 


where  crw  is  the  undisturbed  wall  shear  stress,  while  Ph(<?w)  the  measured  hole- 
pressure  depending  on  cr° .  (2.2"')  is  a  good  approximation  for  (2.2")  or  (2.2')  if 
\crw  —  &w\  is  small.  According  to  Malkus  &  Webster  [17]  \aw  —  <7°  |  can  be  minimized 
by  proper  die  design  with  h/w  >  1.  The  full  instrument  simulation  results  presented 
in  [17]  also  show  that  the  difference  by  using  either  aw  or  cr^  in  predicting  N\  is 
rather  small. 

3.2  Formulation  of  Streamwise  Stress  Gradient 

fti  dg i 

The  general  expression  of  Tn  is  given  by  equation  (3.1).  We  need  to  differentiate  I’ll 
with  respect  to  streamline  coordinate  q\.  Reader  is  referred  to  reference  [17]  for  the  details 
of  formulation.  Here  we  just  simply  present  the  final  expression  of  the  streamwise  gradient 


1  dfn  _  1  dp0  1  d2u i 

hi  dq\  hi  dqi  +  ^ h\  dq\ 


(3.20) 


4.  Error  Cancellation  -  A  Fortuitous  Phenomenon 


In  order  to  make  a  distinction  between  the  modified  hole-pressure  prediction  defined 
by  (2.15)  and  the  original  HPBL  prediction,  we  shall  use  Ph  to  denote  the  MHPR  and 
reserve  Ph  for  the  original  HPBL  prediction  given  in  (2.1).  With  new  notations,  (2.15) 
can  be  re-written  as 


where 


pi 

rH 

p 2 

r H 


Ph^P'h  +  Ph 


i'll  —  T*22  1  9Ti2 
2 f12  h2  dq2 


Tn  —  T72  1  dfn 

2 T12  hi  dqi 


(4.1) 

h2dq2  , 

(4.2) 

h2dq2  . 

(4.3) 
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It  is  easy  to  see  that  Pjj  is  the  path  integral  corresponding  to  the  original  HPBL  prediction 
Ph  (but  they  are  not  equal,  we  will  see  it  later)  and  Pfj  is  the  streamwise  stress  gradient 
term  neglected  by  the  HPBL  theory. 

In  reference  [12],  Baird  et  a 1.  reported  their  observation  on  some  partial-cancellation 
of  the  HPBL  integral  (2.1)  itself.  They  suggested  that  the  success  of  HPBL  prediction  was 
attributed  to  the  cancellation  of  certain  contributions  to  the  integral  in  (2.1)  until  reaching 
the  (channel)  centerline  at  which  point  one  obtain  just  the  contributions  from  viscometric 
flow  (p.642  [12]).  We  believe  that  their  picture  was  incomplete.  By  studying  the  MHPR 
based  on  the  streamline  coordinate  formulation  and  the  two  key  flaws  of  the  HPBL  theory, 
we  found  a  more  important  error  cancellation  phenomenon  in  the  derivation  of  the  HPBL 
prediction,  i.e.  the  error  caused  by  neglecting  the  stress  gradient  term,  can  be 

canceled  out  by  the  error  introduced  through  the  viscometric  flow  assumption  (A5).  By 
using  the  streamline  coordinate  formulation,  we  can  prove  that  this  error  cancellation  is 
exact  for  the  second-order  fluid.  Furthermore,  the  conclusion  of  error  cancellation  is  also 
generalized  to  include  the  viscoelastic  liquids  described  by  the  viscometric  constitutive 
equation  (2.7). 

4.1  Basic  Results  for  Second-Order  Fluid 

In  order  to  show  the  exact  error  cancellation  in  MHPR  for  the  second-order  fluid,  we 
need  some  primary  results  to  proceed. 

Proposition  4.1  For  incompressible,  steady,  creeping  flows  over  a  transverse  slot,  the 
hole-pressure  difference 

Ph  =  -(T22I,.  -  f22|,; )  =  ^ 

is  exact  for  the  second-order  fluid.  Here  N\  is  interpreted  as  the  generalized  viscometric  part  of 
Tn  —  T22,  i.e.  Ni  =  1/172. 

The  proof  of  Proposition  4.1  has  already  been  given  in  the  formulation  in  section  3. 

Proposition  4.2  Let  p°  be  the  pressure  of  Stokes  solution.  The  following  relation 


1  dp0  _  2pcPui  2Tu  dhi_  1  df12 
hi  dqi  h\  dq\  +  hxh2  dq2  +  h2  dq2 


(4.4) 


holds  on  the  hole  centerline  (£ . 

Proof :  Since  p°  is  the  pressure  of  Stokes  solution,  the  equation  of  motion  gives 


Vp°  =  77V  •  A 


(«) 


where  A  =  2D,  and  the  components  of  V  ■  A  can  be  obtained  by  replacing  T  with  A  in 
equation  of  motion.  The  first  equation  of  (a)  becomes 


1  dp° 
h\  dq\ 


h\h2 


d  2  \  &  ,,  :  \  a  a  &h2 

+  AnW2~  22 Wi. 


(6) 


15 


By  using  the  symmetry  conditions  (2.8),  (b)  can  be  easily  reduced  to  (4.4).  0 
Proposition  4.3  For  second-order  fluid,  we  have 

Ph  =  A?./4  -  P3h 


(4.5) 


where 


„„T  (*  fi.  a2a.  i  dfI2 

H  ^  I  i2 


'q$  •'l  T\2  9q2 

Proof :  Substituting  (3.17)  into  (4.2)  yields 

p],  =  f 

H  Jqi  \  h\  dq\  )  2 T12  dq2 


=  P 


2r)Ti 


12  ^92 

2  7  97  ,  1  9T12 


-C^2 


-/ 


2r?T- 


2r]idq2~'1*  Jq .  “*'*  hf  dg?  f12  0?2 


d?2 


=  (il)2  -  Pj  =  JC./4  -  Pf,  0 


(4.6) 


Remark  :  Note  that  Pf  can  be  considered  as  the  error  term  introduced  by  the 
viscometric  flow  assumption  (A5).  According  to  our  numerical  results,  Pj^  is  in  general 
non-zero.  Therefore  Proposition  4.3  indicates  that  in  general  P^  ^  iVf^/4  and  the  error 
included  in  Pf  term  comes  from  the  viscometric  flow  assumption. 

Proposition  4.4  For  the  second-order  fluid,  we  have 

Ph  =  <„/4  -  P'H  .  (4.7) 


Proof  :  By  using  (3.20)  and  (4.4)  in  (4.3) 


P2  =  r*  (  1  dp° . 

H  7,5  V  *i  09 1  09 1 


Tn-T; 


22 


r 


2Ti2 

1  5Ti2  \  Tn  —  T22 


2T12  5/ii 
hi/i2  d?2  h2  5^2  /  2T] 


h2dq2 
h2dq2 
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9“  n  1  oq2 

^22 
dq2 


dq2  -Pf  =  N*J 4  -  Pf 


where  equation  (2.17)  and  Proposition  4.1  have  been  used.  0 

Remark  :  Keep  in  mind  that  Pf  is  actually  the  error  term  caused  by  neglecting  the 
streamwise  stress-gradient,  Proposition  4.4  tells  us  that  for  second-order  fluid  the 
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MHPR  predicts  the  same  result  as  the  HPBL  theory,  namely,  Ph  =  P h  +  P h  —  N°w/ 4  = 
Ph- 

Based  on  Propositions  4.3  &  4.4  we  can  easily  prove  the  following  proposition  regarding 
the  exact  error  cancellation  in  the  MHPR  for  the  second-order  fluids. 

Proposition  4.5  For  incompressible,  steady,  creeping  flows  of  the  second-order  fluids 
over  a  transverse  slot,  we  have  the  following  relation 

p 2  _  p3 

*  H  ~  *H  ‘ 

Proof :  By  simply  combining  Propositions  4.3  &  4.4.  <0> 

Remark  :  The  result  presented  in  Proposition  4.5  is  interesting  and  important.  It 
concludes  that  for  second-order  fluid  the  error  caused  by  one  key  flaw  in  the  original  HPBL 
derivation  is  exactly  equal  to,  hence  is  exactly  cancelled  out  by,  the  error  introduced 
through  another  key  flaw.  In  order  to  make  this  point  more  explicitly  we  rewrite  (4.1)  in 
the  following  form: 

pH-Pk  +  Ph-NiJi  -Ph  +  Ph  —  Ph  ■ 

S  V  *  ^ 

=0 

This  fortuitous  error  cancellation  phenomenon  has  been  unfortunately  hidden  behind 
the  success  of  the  HPBL  prediction  for  a  long  time,  and  now  it  eventually  becomes  clear 
to  us.  Based  on  Proposition  4.5  the  following  conclusion  can  also  be  drawn,  namely: 

For  incompressible,  steady,  creeping  flows  of  the  second-order 
fluid,  the  MHPR  in  path  integral  form,  (2.15)  is  equivalent  to  the 
re-interpreted  HPBL  prediction  in  stress  integral  form,  (2.1'),  due 
to  the  exact  error  cancellation. 

4.2  Generalized  Results  for  the  Viscometric  Model 

Having  obtained  the  basic  results  for  the  second-order  fluid,  we  now  turn  to  the  fluids 
described  by  the  viscometric  constitutive  equation  (2.7),  and  try  to  generalize  the  error 
cancellation  theorem  to  include  more  fluid  models.  Similar  to  the  case  of  second-order 
fluid,  we  begin  by  presenting  some  generalized  results  for  the  viscometric  model. 

Proposition  4.6  For  the  viscometric  model  described  by  (2.7),  the  physical  stress 
components  on  the  hole  centerline  <£  are 

f12 

Tn 

T22 


=  7(7)7 

=  -p  +  {vi  +  i/2)72  -  Vi 


u  1  d2u  1 
h\  dq\ 


.  2  Ul  d2u  1 
-p  +  1^27  +  V\  ■ 


-  V 1 


h\  dq\ 

V  M  a,?  ) 


(4.9) 


(4.8) 
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where  ult  etc.  is  now  the  generalized  Newtonian  solution  associated  with  77(7). 

Proof :  Simply  inserting  (3.9),  (3.11)  and  (3.15)  into  (2.7)  directly  gives  the  above 
results.  0 

Remark  :  (i)  It  should  be  borne  in  mind  that  77,  17,  1/2  are  generally  functions  of 

7.  When  77  =  constant,  u\  =  2rjT  and  1/2  =  —tjT/2,  (4.9)  reduces  to  the  second-order 
fluid  result  (3.16)  &  (3.17).  (ii)  Although  equation  (2.7)  is  called  “ viscometric ”  model, 
the  stresses  still  may  not  be  viscometric  due  to  the  disturbance  of  the  slot,  as  one  can  see 
from  (4.9).  When  flow  is  not  viscometric,  1/1  is  only  the  viscometric  part  of  the  first  normal 
stress  coefficient. 

Proposition  4.7  For  the  viscometric  model  under  Tanner’s  assumptions  (Tl)  to  (T3), 
the  hole-pressure  difference  can  be  expressed  as 


(4.10) 


Proof :  From  (4.9)  and  (2.10)  we  have 


f22k  =  -p°(o,?2a) 

?22l,»  =  -p°(0,?S)  -  Y  (id)2  +  \Jo  v\i  di 

Noting  />°(0,?2)  =  0  on  <£,  (4.10)  can  be  immediately  obtained.  0 

Remark  :  (4.10)  was  first  given  by  Tanner  in  [20].  As  a  matter  of  fact,  (4.10)  has 

another  equivalent  form. 

Proposition  4.8  For  the  viscometric  model  under  Tanner’s  assumptions  (Tl)  to  (T3), 
the  Tanner’s  result  (4.10)  can  also  be  written  as  the  following  equivalent  form: 

•  0 

Ph  =  +  f  /  •  (4-11) 

Proof :  Integrating  (4.10)  by  parts  yields 


where  V\  =  arj.  On  the  other  hand,  integrating  (4.11)  by  paxts  gives  (4.10).  0 

Next,  we  prove  that  results  (4.10)  and  (4.11)  axe  equivalent  to  the  re-interpreted  HPBL 
prediction  (2.1'). 

Proposition  4.9  For  the  viscometric  model  under  Tanner’s  assumptions  (Tl)  to  (T3), 
the  Tanner’s  hole-pressure  difference  results  (4.10)  and  (4.11)  are  equivalent  to  the  re-interpieted 
HPBL  prediction  (2.1').  when  Nx  in  (2.1')  is  interpreted  as  only  the  viscometric  part  of  fn  -T22. 
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Proof :  Substituting  N\  =  apj2,  o  =  777  into  (2.1') 


74*77) 


which  gives  (4.10).  On  the  other  hand,  by  using  the  following  relations 


*  0  ^0 

\Niw  =  \  j  4a*772)  =  J  +  21771*7)  (d) 

0  *  0 
AT.  ry 

2^da  =  2  Jo  ^  dT/  +  V^d^  ^ 

in  (4.11)  we  can  obtain  (2.1').  0 

Now  it  is  the  time  to  consider  the  MHPR  for  the  viscometric  model.  To  proceed,  it  is 
essential  to  have  the  expression  for  stress  gradient. 

Proposition  4.10  For  the  viscometric  model  under  Tanner’s  assumptions  (Tl)  to  (T3), 
the  streamwise  stress  gradients  on  hole  centerline  £  are 


1  dp  _  1  dp° 
dq\  hi  dqi 

hi  dqi  hi  dqi  h\  dq\ 


(4.12) 


where  p  is  the  modified  pressure  given  by  (2.10)  and  p°  is  the  pressure  for  the  generalized 
Newtonian  flow. 

Proof :  Differentiating  (2.10)  with  respect  to  qi  yields 


1  dp 
hi  dqi 


1  dp° 
hi  dqi 


a  d  /tii  dp°\  1  9 

2hi  dqi  \hi  dqi  )  hi  dqi 


vii2 / 2  +  u2i2  -  - 


Keeping  in  mind  that  uj,  h  1,  7,  dp0 /dqi  are  even  functions  of  qi  on  £,  consequently  the 
last  two  terms  in  above  equation  are  all  zero.  Then  the  derivation  of  is  similar  as 

in  the  second-order  fluid.  <> 

Remarks  :  Comparison  of  (4.12)  with  (3.20)  tells  us  that  the  streamwise  stress 
gradient  of  the  viscomertic  model  is  the  same  as  that  of  the  second-order  fluid.  Therefore, 
the  proof  of  the  error  cancellation  theorem  for  the  viscometric  model  will  be  very  similar  to 
that  of  second-order  fluid.  Since  we  have  assumed  that  the  velocity  field  for  the  viscometric 
model  (2.7)  are  the  same  as  the  generalized  Newtonian,  flow  field,  Proposition  4.2  is  still 
useful,  while  Propositions  4.3  &  4.4  need  to  be  slightly  modified  as  follows. 


Proposition  4.11  For  the  viscometric  model  under  Tanner’s  assumptions  (Tl)  to  (T3), 
we  have  the  following  results  for  the  path  integral  of  MHPR 


where 


1  Pw 

P'h  =  <„/2  -  5  J  vxiH  ~  P3„ 

P2h  =  KJ2  ~  \  £"  "i  Hi  -  P'h 
Ph  —  Pr 


3  _  fq 2  Ui  d2U1  V\  dfi2  , 

H  Jai  h\  d<ll  Ti2  92 


(4.13a) 

(4.136) 

(4.13c) 


(4.6a) 


The  proof  of  Proposition  4.11  is  straightforward  and  has  been  omitted  here.  On  the 
basis  of  Proposition  4.11,  one  can  easily  prove  the  following  result  regarding  the  equivalent 
expressions  for  Ph- 

Proposition  4.12  For  incompressible,  steady,  creeping  flows  of  the  viscometric  model 
(2.7)  with  Tanner’s  assumptions  (T1)-(T3),  there  are  four  equivalent  expressions  that  predict  the 
hole-pressure  error,  Ph-  They  are  the  MHPR  (2.15),  the  Tanner’s  results  (4.10)  &  (4.11),  and 
the  re-interpreted  HPBL  formula  (2.1'). 

Remarks  :  From  section  2  we  know  that  the  MHPR  is  derived  from  the  equations  of 
motion,  without  any  assumptions  on  the  fluid  model.  Therefore  (2.15)  is  model  indepen¬ 
dent  and  hence  the  more  general  expression  for  the  hole-pressure  error.  For  the  viscometric 
model,  all  the  other  three  relations  can  be  derived  from  the  MHPR  under  the  Tanner’s 
assumptions  (Tl)-(T3).  It  is  also  easy  to  see  that  when  77 ,  tq  =constant,  all  the  four 
expressions  give  the  same  result  Ph  =  Ar1°u,/4. 


5.  Conclusions 

The  correct  application  of  the  HPBL  theory  in  the  hole-pressure  problem  involves 
the  modified  hole-pressure  relation  (MHPR)  in  path  integral  form.  Some  important 
mathematical  properties  of  the  MHPR,  such  as  the  singularities  and  the  existence  of  the 
path  integral,  were  studied. 

To  investigate  the  two  key  flaws  of  the  HPBL  theory,  it  is  sufficient  to  consider  the 
second-order  fluid.  By  studying  the  MHPR  in  terms  of  the  streamline  coordinate  for¬ 
mulation,  we  found  the  fortuitous  phenomenon  of  error  cancellation  in  the  derivation  of 
HPBL  formula.  For  second-order  fluid  and  Tanner’s  viscometric  model  under  certain  as¬ 
sumptions,  we  have  proved  analytically  that  the  error  cancellation  is  exact  and  the  MHPR 
is  equivalent  to  the  re-interpreted  HPBL  prediction.  We  also  proved  that  there  axe  four 
equivalent  expressions  that  predict  the  same  hole-pressure  error  Ph  for  Tanner’s  visco¬ 
metric  model.  The  theoretical  results  of  this  paper  provided  a  complete  explanation  for 
the  fortunate  success  of  the  HPBL  prediction  and  the  controversy  in  the  derivation  of  the 
HPBL  formula. 
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